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The influence of a magnetic field on superconductivity is usually described either 
phenomenologically, using Ginzburg-Landau theory, or semiclassically using Gor'kov 
theory. In this article we discuss the influence of magnetic fields on the mean-field 
theory of the superconducting instability from a completely quantum mechanical 
point of view. The suppression of superconductivity by an external magnetic field is 
seen in this more physically direct picture to be due to the impossibility, in quan- 
tum mechanics, of precisely specifying both the center-of-mass state of a pair and 
the individual electron kinetic energies. We also discuss the possibility of novel as- 
pects of superconductivity at extremely strong magnetic fields where recent work has 
shown that the transition temperature may be enhanced rather than suppressed by 
a magnetic field and where a quantum treatment is essential. 



1 



I. INTRODUCTION 



The relationship between superconductivity and magnetic fields is both of practical im- 
portance in the design of superconducting devices and of fundamental importance to the 
superconductivity phenomenon. In the absence of an external magnetic field superconduc- 
tivity is associated with the pairing of time-reversed electron states. As we discuss in detail 
below, magnetic fields break time-reversal-invariance symmetry and frustrate this pairing. 
For sufficiently weak external magnetic fields superconductors prefer to completely expel any 
external magnetic fiux (the Meissner effect) in order to avoid this frustration. At stronger 
magnetic fields type-II superconductors, which are used in the construction of supercon- 
ducting magnets, can form a mixed state in which superconductivity coexists with magnetic 
fiux. Superconductivity in the mixed state is usually described in terms of Ginzburg-Landau 
theory p6[ which predicts a decrease in the temperature to which superconductivity can sur- 
vive (Tc) proportional to the external magnetic field strength. For sufficiently weak external 
fields and temperatures close to Tc Ginzburg-Landau theory was derived microscopically by 



Gor'kov [|T2|. Gor'kov's theory has a wider range of validity than Ginzburg-Landau theory 
and predicts ||13|Jl^ ,p| that Tc decreases monotonically with increasing magnetic field and 
is eventually driven to zero. However Gor'kov's theory treats the magnetic field in a semi- 
classical approximation which is not valid when the temperature is sufficiently low and the 
disorder is sufficiently weak that the Landau quantization of motion in planes perpendicular 
to the field direction becomes important. In the past few years , following seminal work by 



Rasolt, Tesanovic and collaborators, it has been realized that (at least within the stan- 
dard mean-field-theory known to be accurate at weak magnetic fields) superconductivity can 
survive to arbitrarily strong magnetic fields once Landau quantization is accounted for. In 
this article we discuss the superconducting instability in a magnetic field from a completely 
quantum-mechanical point of view. We explain how the results of Ginzburg-Landau theory 
and Gor'kov theory can be understood in terms of the microscopic quantum mechanics of 
charged particles in a magnetic field and why Gor'kov theory can fail at sufficiently strong 
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fields. 



II. Tc AT ZERO MAGNETIC FIELD 

It is useful to begin by discussing the familiar implicit equation for in the absence of 



a magnetic field : 



(Energies are measured from the chemical potential u and fl is the volume of three- 
dimensional systems or the area of two dimensional systems. The Fermi energy, ep = mVp/2 
is the zero temperature limit of fi.) This equation is for the usual BCS model with attractive 
interactions of constant strength V. All the discussion in this article will be in terms of this 



simple model . The prime on the sum over wavevectors denotes the usual separable high 
energy cutoff requiring both electron energies to be within of the Fermi level. The nu- 
merator of the factor in square brackets in Eq.( expresses through the Fermi occupation 
numbers the requirement that the pairing come either from electrons outside the Fermi sea, 
as in the Cooper problem, or from holes inside the Fermi sea. Note that this factor vanishes 
at finite temperature, and even at T = for P 7^ 0, when + ep is near zero. In a su- 
perconductor a bound state occurs for the relative motion of electrons in a Cooper pair and 
the temperature at which the bound state first occurs, Tc, depends on the center-of-mass 
(COM) momentum of the pair, P, as we discuss in the following paragraph. 
Defining an effective pairing density-of-states by 

^"(e --P^T)^^ ±[1 - f{e^) - f{ep_m^ - 6,^ - ep_^\ (2.2) 

k 

the Tc equation can be rewritten in the form 

l = V z/P(e : P,T)/e. (2.3) 



At P = and T = z/P(e) = e{2E+ - |e|)(e/|e|)z/(2e)/2 where z/(e) is the single-electron 
density of states per spin. At finite P and T z/^(e : P, T) is reduced toward zero for 



|e| <~ E~ = sup{kBT, VpP) because of the combination of Fermi factors appearing in but 
otherwise is nearly constant. (Here Vp is the Fermi velocity. Low energy pairs tend to be 
composed of states on opposite sides of the Fermi energy for P ^ 0. See Fig.( |I]). The 
pairing densities-of-states for two-dimensions (2D) and three-dimensions (3D) are reduced 
as a consequence when |(e — 2ei?)/VFP| < 1.) For T = the reduction in z/^ is illustrated in 
Fig.( I). For |(e - 2eF)/VFP\ < 1 

for 3D and 

uP{e : P) = ^(1 - - cos-\{e - 2e^)/r^P)). (2.5) 

for 2D. For weak coupling {E~ <^ E^) the Tc equation reduces to E~ ~ £'+exp(— 1/A) 
which will have no solution once VpP exceeds ~ kBTc{P = 0). (A = Vu^O). It is assumed 
that z/(e) is nearly constant over the energy range E~^.) Later we will relate this result for 
the dependence of Tc on the COM momentum of the Cooper pair directly to the dependence 
of Tc on an external magnetic field. 



III. PAIR STATES IN A MAGNETIC FIELD 

Note that the equation depends both on the COM state of the pair and, through the 
Pauli-exclusion-principle requirements expressed by Fermi factors, on the states of the indi- 
vidual electrons making up the pair. The states of a pair of electrons may be described either 
in terms of COM and relative motion states or in terms of the individual electron states. In 
the absence of a magnetic field this connection is trivial. To describe superconductivity in a 
magnetic field quantum mechanically we must start by discussing the relationship between 
these two descriptions in a magnetic field. The Hamiltonian |^ for two non-interacting 
electrons , h, is 

h = -^(-*Wi + -A{n)f + ^{-thV2 + -A{r2)f (3.1) 
2m c 2m c 
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or 

Here we have assumed a gauge where the vector potential is hnear in the coordinates, 
M — 2m , 11 — m/2, R = {fi + /2 and f — fi — r2. Notice that the charge appearing 
in the center of mass term is 2e while the charge appearing in the relative motion term is 
e/2 so that both relative and center of mass kinetic energies (KE's) are quantized in the 
same units as for individual electrons, hujc — eB/mc. (The individual electron eigenvalues 
measured from the chemical potential are ejv = hu!c{N + 1/2) — fi = hu!c{N — Np). Np 
is the Landau level index at the Fermi level.) In the Landau gauge the eigenfunctions for 
individual electrons are well known: 

^N,x{fi} = exp{-iXyi/e^)Mi^i - ^))/V^ (3-3) 

where Ly is the length of the system in the y direction, i = (he/ eB)^/'^ is the magnetic length, 
and 4>n{x) is a one-dimensional harmonic oscillator eigenstate for mass m* and frequency 
cUc- The expressions for the center-of-mass and relative eigenstates, and are identical 
except that the characteristic lengths are scaled to account for the changes of charge and 
mass. (The effective magnetic lengths are = i/\/2 and = \/2i for the center-of-mass 
and relative eigenstates respectively.) 

In the lowest Landau level (j)N=o{x) ~ exp(— x^/4£^) so that 

tpo,x+Y/2{ri)ipo,x-Y/2{r2) = '0o,x(^)V'S,y(O- (3-4) 

The relationship is easily generalized to higher Landau levels by writing the Hamiltonian in 
terms of ladder operators. 

h = hu!c{a\ai + 0202 -f- 1) = hurdci'j^ciR + alar + 1) (3.5) 

and noting that = (ci + a2)/V^, and = (oi — a2)/\/2- Here aj = -^^{T^xj ~ iT^yj)-, ^-^id 
TTj- = —ihVj -\- ^Aj. It follows that 
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N+M 

i^N,x+Y/2{ri)ijM,x-Y/2{r2) = ^ Bf'^'^f^xiRWN+M-jA"^ (3-6) 

j=0 

where 

^ J\{N + M - j)\N\M\ „2 ^ i-r-^ 

^ ^ 2^+^^^ ^ ^0 (j -m)!(Ar + m-j)!(M-m)!m! ^ ' ^ 

Note that both left and right hand sides of Eq.( |3]^) are manifestly eigenstates of h with 
eigenvalue huJc{N + M + 1). The coefficients B^'^^ give the amplitude for having KE hujc{j + 
1/2) in the center-of-mass motion ( and huJc{N + M — j + 1/2) in the relative motion) when 
the individual particles have definite KE's huJc{N + 1/2) and huJc{M + 1/2). 

The coefficients appearing in the unitary transformation between the two sets of two- 
particle eigenstates, {B^'^^} will play a central role in the discussion below. Note that the 
transformation is block-diagonal with no mixing between eigenstates of different total kinetic 
energy. The completeness of either set of eigenstates implies the following identities: 

K 



E B^^-'^Bf'^-'' = (3.8) 

N=0 



K 

Z^Bj Bj =On',n (3.9) 

j=0 

Since the center of mass kinetic energy does not commute with the individual particle ki- 
netic energies the center of mass is necessarily uncertain if the individual particle states 
are known precisely. Conversely, for given center-of-mass and relative state kinetic energies 
the individual particle kinetic energies are necessarily uncertain. Given j and the relative 
motion eigenstate, or equivalently j and the total kinetic energy index K, gives 
the normalized probability distribution for the individual electron states with the same total 
kinetic energy. Explicit expressions for small j are easily obtained from Eq.( |3.7] ): 

\Br-r-U'^\ (3.10) 



I TjN,K-Ni2 
l-Dl I 



'^^ ' =2^\n) 2K{K-1) ■ ^^-^^^ 
For K » \k\, {k = N - M) it can be shown |[ that 

j^iK+k)/2,iK-k)/2 ^ ^l_y/^^±_f^H,{k/V2K) expi-e/4K) (3.13) 
where Hj is a Hermite polynomial. 



IV. Tc IN A MAGNETIC FIELD 



The implicit equation in a magnetic field [10| is completely analogous H to the B = 



equation (Eq.( |2ri|) ) cited at the beginning of this article: 



As in the B = case Tc depends on the Cooper pair state. As we discussed previously 
the superconducting Tc decreases with \P\ for B = 0. For B ^ 0, Tc is independent of 
the guiding center quantum number X for the Cooper pair. The fact that instabilities 
occur simultaneously in a macroscopic number of channels is responsible for the dimensional 
reduction which causes superconducting fiuctuations to be qualitatively altered by a 
magnetic field. The superconducting instability still depends, however, on the Landau level 
index of the Cooper pair. We first examine the weak field limit where ksT ^ Tiuj^. In this 
limit the sums over Landau levels may be replaced by integrals and Eq.( [4.1| ) becomes 

(We've noted that z/(0) = l/(27r£2^u;c) and ir+ is the maximum kinetic energy index allowed 
by the high energy cutoff.) To understand why superconductivity is suppressed by weak 
magnetic fields it is sufficient to consider the T = limit. The Landau levels with indices 
[K + k)/2 and {K — k)/2 are on the same side of the Fermi level and can contribute to the 
pairing only if \k\ < \K — 2Nf\- (See Fig.( ^).) For a given center-of-mass index j of the 
Cooper pair and a given total kinetic energy the probability of finding both members of a 



Cooper pair on the same side of the Fermi energy {ep = fJ'{T = 0)) , is necessarily less than 
one. In Fig.( ^ we plot 

Pj{K) = ^[1 - e{2NF -K-k)- e{2NF -K + (4.3) 



for j = and Np = 12.5 against K. From Eq.( |3.13D we see that most of the contribution 



to Pj{K) comes from \k\ <~ y (j + 1/2)K. The logarithmic divergence of the integral over 
K in Eq.( [4.2|) which guarantees a solution is therefore cutoff since Pj{K) will fall to zero 



for \K — 2Np\ <~ Y (2j + 1)Nf- It follows that solutions at T = exist only if 

hcoc <~ {kBT,f/{2j + l)eF. (4.4) 

The superconducting instability is suppressed most weakly for Cooper pairs with j = 0, i.e. 
for COM in the lowest Landau level, in agreement with Ginzburg-Landau theory. 

At zero magnetic field the superconducting instability occurs first for COM momentum 
P = 0; the pairing of time-reversed single-particle states guarantees that all pairs are allowed 
by the Pauli exclusion principle at T = even if their energies are very close to the Fermi 
energy. In a magnetic field time-reversal-symmetry is broken so that time reversed pairs 
of single particle states no longer exist. The kinetic energy eigenstates in a magnetic field 
are usefully thought of as having a definite magnitude of momentum corresponding to the 
quantized kinetic energy but completely uncertain direction of momentum since they are 
executing circular orbits. For definite COM and relative kinetic energies, en and e^, the 
mean-square difference in individual electrons kinetic energies is 

((ei - e2f)9 = 2eRer (4.5) 

The average here is over the angle between the COM and relative momenta which is 
completely uncertain in a magnetic field. This classical root-mean-square energy differ- 
ence agrees with the energy width of the quantum mechanical distribution function dis- 
cussed above. When the mean energy of the pair is within this width of the Fermi en- 
ergy contributions to pair formation are suppressed by the Pauli exclusion principle. For 
= hu!c{j + l/2)^er-~2ei? the resulting low-energy cutoff is 
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E- =~ 2^nuJc{j + l/2)eF = VpPj. (4.6) 

In Eq.( ^) hPf/Am = hujc{j + 1/2) so that Pj is the 'quantized' magnitude of the COM 
momentum. We see from this discussion that pairing in COM Landau level j in a magnetic 
field is very similar to pairing at COM momentum Pj in the absence of a magnetic field. 

The above discussion explains from a quantum mechanical point of view the familiar 
suppression of superconductivity by a magnetic field in the weak field regime where the 
discretization of allowed kinetic energies transverse to the magnetic field is washed out 
either by temperature of disorder. In clean 2D systems the Landau level structure becomes 
important in the thermally averaged density of states for huc >~ ksT; in 3D systems the 
free motion along the magnetic field partly obscures the Landau level structure and the 
strong field regime is reached only for huc >~ ^/NphsT . In the strong field regime the 
density-of-states has strong peaks and the chemical potential tends to be pinned to these 
peaks. It is these density of states peaks which can reverse the decrease of T^. with field 
and lead to a peculiar regime where increases with field. As the strong field limit is 
approached the Landau level at the Fermi energy contributes more strongly to the sum 
in Eq.( [4.11 ). One immediate effect apparent even at comparatively weak fields ^ is the 
decrease in Tc for j odd. ( For COM j odd the probability of pairs occupying the same 
Landau level is zero.) Magnetooscillations [^,^,0 in T^, and in all properties of the mixed 



state 1^ of the superconductor occur as Landau levels pass through the Fermi level. These 
oscillations have been observed experimentally |^,17| and are not yet understood in complete 
detail. 

At extremely strong fields a regime can be reached where only electrons in the Landau 
level at the Fermi energy contribute importantly to the pairing. In this limit (for 2D systems) 



Tc reaches a maximum when the Landau level is half full ||2T| and Eq( |4.1|) reduces to 



T., = ^|sf-^-|2 (4.7) 

Note that Tcj is proportional to magnetic field strength. In the extreme quantum limit all 
electrons are in the lowest Landau level and Np = 0. Since the maximum value of j is 



2Nf it happens that superconductivity occurs in the j = channel just as in the weak 
magnetic field limit. This similarity in the nature of the superconducting order in the weak 
and infinitely strong field regimes suggests that no novel behavior can occur at intermediate 
fields. The suggestion is misleading as we can see by looking at the case where Np 7^ 0. The 
maximum COM kinetic energy channel for the Cooper pair is 2Nf and pairing can occur in 
any even j channel. From the expression for B^^'^^ we find that in this case Tc tends to be 
larger for j close to either its minimum or maximum values and is always equal for j = 
and j = 2Nf- (See Table I) This result can be understood by calculating the probability 
that two electrons of the same energy ep but with completely uncertain relative orientation 
of momentum will have a given COM kinetic energy, e^. Averaging over angles it is easy to 
show that 

P{eR) = -{en{2eF - en))-''^ (4.8) 

TT 

which is peaked near the minimum and maximum possible values for e/j. Thus in the 
extremely strong field regime there is the possibility of unusual superconducting states in 
which Cooper pairs are in states with elevated kinetic energies. In mean-field theory the 
vortex-lattice state is found ||2T|] to have j > and to have associated unusual properties 



including the possibility of having several vortices per period of the lattice. 



V. CONCLUDING REMARKS 

In this article we have discussed how the suppression of superconductivity by a magnetic 
field can be understood completely microscopically in terms of the quantum mechanics of 
pairs of particles in a magnetic field. The results obtained in this way are equivalent to those 
obtained by Ginzburg-Landau theory and Gor'kov theory in their ranges of validity. The 
suppression is related to the quantum uncertainty in the kinetic energies of the individual 
electrons making up a Cooper pair of definite COM kinetic energy. We have also discussed 
how the suppression can be overcome by the enhancement of the density of states near 
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the Fermi energy which occurs for sufficiently strong magnetic fields in clean samples and 
explained why the Cooper pair wavefunction can be unusual in this regime. Ginzburg- 
Landau theory is not valid in the regime of strong-field superconductivity except for the 
case where pairing occurs entirely within the = Landau level. 

We have restricted our attention here to aspects which follow directly from the quantum 
mechanics of pairs of charged particles in a magnetic field and the reader should be aware that 
many other issues arise, some parasitically, especially when considering superconductivity 
in extremely strong magnetic fields. For example, in our discussion we have, for the sake of 
definiteness, taken the electron gf— factor to be zero; a non-zero g-factor will affect results 
at strong fields p^. For the sake of our discussion here we have also assumed that the 
standard mean-field theory of superconductivity which leads to the expressions for we 
have employed and which is known to be relaible at weak fields can also be used in the 
strong field regime. It is certain that this is not entirely correct especially in the case of 
2D case although we beheve that the considerations discussed here are still essential for 
the physics in that regime. 
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Superconductivity Consortium through D.O.E. grant no. DE-FG-02-90ER45427 and in part 
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FIGURES 

FIG. 1. Pauli blocking of low-energy pair states at finite center-of-mass momentum, P, of the 
pair. In terms of the relative momentum of the pair the Fermi surfaces for the two electrons 
composing the pair are displaced by P. The pair must be composed of unoccupied electron states 
or hole states. At zero temperature the allowed values of relative momentum are either inside both 
Fermi surfaces or outside both Fermi surfaces. The shaded regions where the energies are close to 
the Fermi energy are forbidden. 

FIG. 2. Pairing density of states at finite COM momentum for three-dimensions (solid line) 
and two-dimensions (dashed line). 

FIG. 3. Probability of having individual electron kinetic energies fiLOc{K + k + l)/2 and 
huJc{K — k + l)/2 given center-of-mass kinetic energies hujc{j + 1/2) and total kinetic energies 
huJc{K + 1/2). The probabilities are represented by the vertical lines at even integer values of k. 
(k must be even when K is even and odd when K is odd.) The results shown here are for K = 30 
and j = 0. For Np = 12.5, i.e. for the first 12 Landau levels occupied, the two single-particle 
states are both occupied or both empty only for k = 0, k = ±2 and k = ±4. Larger values of k, 
for which the probability is indicated by a dashed line, are Pauli blocked and cannot contribute to 
pairing in a j = center-of-mass state. For this case the probability that the two single-particle 
states will be on the same side of the Fermi energy is Po{K = 30) = 0.6384. The solid line which 
envelopes the probabilities is the large K expression Eq.( p.l3|) . 

FIG. 4. Pj against K for j = and Np = 12.5, i.e. for the first twelve Landau levels occupied. 
For the case K = 30 Pj is given by the sum of the probabilities indicated by the solid lines in 
Fig.( I). 
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TABLES 

TABLE 1. IB^'^'^^I'^ for Np = 0, 1, 2, 3. 



J = J = 2 j = 4 j = 6 

10 

1/2 1/2 

3/8 1/4 3/8 

15/48 3/16 3/16 15/48 
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